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Liquid rocket engines using oxygen and methane propellants are being considered by the 
National Aeronautics and Space Administration (NASA) for future in-space vehicles. This 
propellant combination has not been previously used in flight-qualified engine systems 
developed by NASA, so limited test data and analysis results are available at this stage of 
early development. As part of activities for the Propulsion and Cryogenic Advanced 
Development (PC AD) project funded under the Exploration Technology Development 
Program, the NASA Marshall Space Flight Center (MSFC) has been evaluating capability to 
model combustion performance and stability for oxygen and methane propellants. This 
activity has been proceeding for about two years and this paper is a summary of results to 
date. Hot-fire test results of oxygen/methane propellant rocket engine combustion devices 
for the modeling investigations have come from several sources, including multi-element 
injector tests with gaseous methane from the 1980s, single element tests with gaseous 
methane funded through the Constellation University Institutes Program, and multi-element 
injector tests with both gaseous and liquid methane conducted at the NASA MSFC funded 
by PC AD. For the latter, test results of both impinging and coaxial element injectors using 
liquid oxygen and liquid methane propellants are included. Configurations were modeled 
with two one-dimensional liquid rocket combustion analysis codes, the Rocket Combustor 
Interactive Design and Analysis code and the Coaxial Injector Combustion Model. Special 
effort was focused on how these codes can be used to model combustion and performance 
with oxygen/methane propellants a priori , and what anchoring or calibrating features need 
to be applied, improved or developed in the future. Low frequency combustion instability 
(chug) occurred, with frequencies ranging from 150 to 250 Hz, with several multi-element 
injectors with liquid/liquid propellants, and was modeled using techniques from Wenzel and 
Szuch. High-frequency combustion instability also occurred at the first tangential (IT) 
mode, at about 4500 Hz, with several multi-element injectors with liquid/liquid propellants. 
Analyses of the transverse mode instability were conducted by evaluating injector 
resonances and empirical methods developed by Hewitt. 


Nomenclature 

= droplet size constant used in CICM 
= atomization constant used in CICM 
= characteristic exhaust velocity, ft/sec 
= diameter, inches 

= mean droplet diameter used in CICM, microns 
= Rupe mixing efficiency 
= frequency, Hz 
= inner diameter, inches 
= length, inches 

= geometrical chamber length (injector face to geometric throat), inches 


* Engineering Specialist, ER32, Mail Stop ER32, NASA MSFC, Senior Member AIAA. 
1 Aerospace Engineer, ER32, Mail Stop ER32, NASA MSFC, Member AIAA. 
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characteristic chamber length, inches 
jet stripping rate used in CICM 
oxidizer-to-fuel mixture ratio 
pressure, psia 

mass median liquid droplet diameter, microns 

relative velocity, ft/sec 

axial location, ft 

velocity, ft/sec 

efficiency 

absolute viscosity, lbm/ft-sec 
density, lbm/ft 2 3 
surface tension, lbf/ft 


Subscripts 

b 

c 

c* 

g 

ie 

Isp 

j 

L 

m 

mix 

vap 


barrel (cylindrical section) 
chamber 

characteristic exhaust velocity 
gas 

injector end 
specific impulse 
jet 

liquid 

measured 

mixing 

vaporization 


I. Introduction 

N ASA has been considering use of oxygen and methane propellants for propulsion of future space exploration 
vehicles for several years. 1,2 Unfortunately, engine systems with oxygen and methane propellants have never 
been previously certified or flight-qualified by NASA, so test data and analysis results are limited at this stage of 
preparation. NASA has funded two efforts in the past 40 years for thrust chamber development with 
oxygen/methane propellants. During the late 1970s through the 1980s, NASA examined oxygen/methane 
propellants for high-chamber pressure boost phase engine applications. Methane as a booster engine propellant was 
considered a compromise between the density of a hydrocarbon like kerosene and the performance of hydrogen. 
Thrust chamber development activities included the design, development and test of shear and swirl coaxial element 
injectors. 3 " 7 Since the methane was to be used as combustion chamber and nozzle coolant, it was delivered to the 
injectors in gaseous form. High performance (>98% rj c *) was achieved, although high-frequency combustion 
instability was observed with some configurations. 7 The liquid oxygen (L0 2 ), gaseous methane (gCH 4 ) propellant 
conditions were recently revisited in a hot-fire demonstration at the NASA Marshall Space Flight Center (MSFC), 
with similar performance results. 8 Also, as part of university-level research investigations funded by the NASA 
Constellation University Institutes Program (CUIP), a single-element thrust chamber was tested with shear and swirl 
coax element injectors, which included measurement of local heat flux. 9 

More recently, NASA has considered oxygen/methane propellants for in-space rather than booster applications. 
Spacecraft propulsion systems using oxygen/methane propellants have several advantages over systems that use 
hypergolic propellants, including reduction in spacecraft mass due to higher engine specific impulse, the ability to 
use common oxygen tanks for other subsystems (such as life support), and the possibility to produce these 
propellants in situ , thus reducing the need to carry them to extraterrestrial surfaces. 1 In addition, ground crews and 
astronauts can operate in inherently safer environments when manipulating oxygen and methane propellant systems 
compared to hypergolic propellant systems. 1 For these low to moderate chamber pressure engines, where 
turbopump-cycles are less likely to be used, the methane is less likely to be used as thrust chamber coolant, and thus 
would be delivered to the injector in cryogenic liquid form (LCH 4 ). 1,10 Since design, analysis, and experimental 
studies with LCH 4 to date are even less extensive than with gCH 4 , NASA has recently funded, through the 
Propulsion and Cryogenic Advanced Development (PC AD) project under the Exploration Technology Development 
Program at the NASA Glenn Research Center (GRC), several programs for L0 2 /LCH 4 propellant thrust chamber 
development. Four separate commercial thrust chamber programs have been or are being conducted, 10 ' 14 while 
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NASA has conducted additional in-house development activities at the MSFC. 15-18 Performance with LCH 4 
injectors has generally been moderate (rj c *^9 3-96%) with instances of excessively noisy combustion, low frequency 
combustion instability (chug), and high-frequency transverse mode combustion instability. 

As part of these activities, the NASA MSFC has also been funded by PC AD to evaluate the liquid propellant 
rocket engine analysis tools for combustion of oxygen and methane propellants, especially combustor performance 
and combustion stability. The outcome of these evaluations will be an assessment of the readiness to conduct future 
combustion analyses with these propellants, along with the definition of anchoring techniques and recommendations 
of code improvements where applicable. The current one-dimensional combustion and performance codes and the 
combustion Computational Fluid Dynamic (CFD) analysis codes were originally developed largely for other 
propellants and have never been adequately validated for use with oxygen and methane propellants. 

The performance evaluation task at NASA MSFC was divided into 3 subtasks: 1) selection of appropriate 

modeling data sets, 2) evaluation of one-dimensional combustion and performance analysis codes, and 3) evaluation 
of a combustion CFD code. The first subtask gathered existing methane data and selected appropriate data sets that 
can be used to most effectively evaluate the analysis tools. In the second subtask, several of these configurations 
were modeled using the one-dimensional analysis codes Rocket Combustor Interactive Design and Analysis 
(ROCCID) 19 and the Coaxial Injector Combustion Model (CICM). 20 For most of the selected cases, anchoring or 
calibrating values of inputs in the analyses codes were varied to evaluate sensitivities of the analysis techniques and 
to determine the best anchoring values for future analyses. In the third subtask, several of the single-element 
chamber configurations have been and are being analyzed with Loci-CHEM, a CFD code under development at the 
NASA MSFC. 21 

The combustion stability modeling task was a more recent addition due to the growing realization of potential 
combustion instability issues with oxygen and methane propellants, especially with liquid methane. The effort has 
focused to date on evaluation and explanation of the NASA MSFC in-house test results, some of which is reported 
in this paper, but is planned to expand to the contractor programs in the future. 

This paper describes the results of the combustion stability and one-dimensional performance analyses 
conducted to date on several hardware configurations tested with gaseous and liquid methane. Generally, 
descriptions of the hardware are not provided and the reader is encouraged to examine the cited references for this 
information. Also, combustion CFD results are not reported in this paper, but results of ongoing CFD analyses are 
expected to be reported in future papers. 

II. Liquid Oxygen and Gaseous Methane 

Liquid oxygen and gaseous methane combustors were investigated by NASA in the past and have undergone 
more development than combustors with liquid methane. The use of gaseous methane in the injector provides the 
potential for high performance, while combustion instability has occurred in only one instance. Because the 
propellant conditions are liquid and gas, the coaxial element geometry, with the liquid injected in either shear or 
swirl form, is preferred, which brings to bear decades of development from oxygen/hydrogen injector systems. 

A. Performance 

Tests of single-element injectors with shear and swirl coaxial elements and liquid oxygen/gaseous methane 
propellants were conducted at the Pennsylvania State University funded by the NASA CUIP. 9 The original purpose 
of these tests was to measure profiles of circumferential and axial local heat flux, similar to results gathered for 
coaxial elements with liquid oxygen/gaseous hydrogen propellants. ’ ’ These heat flux data may also be used to 
examine the axial profile of energy release rate, since the energy release in the single element injector is likely to be 
similar to the heat flux profile. The overall tj c * can also be measured and compared to calculations. Injector 
geometry and test data for the single-element hardware are summarized in Table I. For all these tests, D c was 1.0”, 
L c was 14.49” with L b of 14.00”, and contraction ratio was 5.38. 

A comparison of performance for various element types is shown in Fig. 1. Note that tj c * for the single element 
experiments were not reported in Ref. [9], but have been calculated for the current analyses. The ?j c * shown in Fig. 1 
do not include heat loss to the chamber, which can provide a substantial increase to the efficiency, on the order of 2- 
4%. Thus, the measured c* without heat loss is likely incorrect since the performance is already approaching 100%. 
However, since all tests were run in the same facility, in close succession and with the same measurement system, 
the values shown in Fig. 1 are believed to be representative and reasonable indicators of the differences between the 
element types. 
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Table I. Single-element geometry and test data from NASA CUIP task at The Pennsylvania State University. 9 


Case Number 


9 

10 

11 

1 

12 

2 

3 

4 

Coaxial Element Type 


Shear 

Shear 

Shear 

OXIDIZER 

Shear 

Shear 

Shear 

Shear 

Swirl 

post characteristics 

__ 

Recessed 

Recessed 

Recessed Recessed 

Recessed 

Rush 

Recessed 

Rush 

metering orifice ID {1} 

in. 

0.092 

0.092 

0.092 

0.092 

0.092 

0.117 

0.117 

3 x 0.020 x 0.332 

post exit ID 

in. 

0.225 

0.225 

0.225 

0.225 

0.225 

0.156 

0.156 

0.244 

exit taper angle 

deg. 

5.8 

5.8 

5.8 

5.8 

5.8 

0 

0 

0 

post exit thickness 

in. 

0.015 

0.015 

0.015 

0.015 

0.015 

0.020 

0.020 

0.020 

post length 

in. 

3.775 

3.775 

3.775 

3.775 

3.775 

1.875 

1.875 

2.225 

post recess 

in. 

0.225 

0.225 

0.225 

FUEL 

0.225 

0.225 

0 

0.156 

0 

annulus length 

in. 

0.500 

0.500 

0.500 

0.500 

0.500 

0.275 

0.275 

0.100 

gap width 

in. 

0.024 

0.024 

0.024 

TEST DATA 

0.024 

0.024 

0.029 

0.029 

0.020 

Test Number 

- 

15 

23 

27 

20 

14 

11 

51 

42 

chamber pressure 

psia 

291 

595 

795 

1006 

1199 

984 

1006 

1016 

mixture ratio 


2.84 

2.98 

3.00 

3.00 

2.99 

3.09 

3.00 

3.00 

fuel manifold temperature 

R 

510 

515 

515 

516 

509 

515 

512 

518 

oxidizer flowrate/e lenient 

Ibm/sec 

0.17 

0.35 

0.47 

0.59 

0.70 

0.60 

0.59 

0.58 

TIC* {2} 

% 

(94.8) 

(96.4) 

(96.6) 

(97.1) 

(96.2) 

(94.1) 

(97.0) 

(98.4) 


Notes: 

{1) swirl coaxial metering is # slots x slot width x slot height 

{2) PSU CUIP performance are estimates without heat loss to combustion chamber 
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Figure 1. Relative performance of single-element injectors from Ref. [9]. Note heat loss is not included. 
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Configurations listed in Table I were modeled using the one-dimensional analysis codes ROCCID 19 and CICM. 20 
The ROCCID code was capable of analyzing both elements. CICM is currently capable of analyzing only shear 
coaxial injection elements with liquid propellant in the center jet and gaseous propellant in the annulus. 

Initial modeling iterations were conducted using pre-defmed calibrating parameters. For the ROCCID code, 
such critical parameters include atomization lengths and droplet sizes, which can be selected from a list of several 
model-provided analytical or empirical calculations (or defined directly by the user or modified in the CONTROL 
file). 19 For CICM, four critical parameters are available for calibration: the atomization coefficient C A in the cup 
and the chamber for calculation of the jet stripping rate M A , shown in Eqn. (1), and the droplet size production 
coefficient B A in the cup and the chamber for the calculation of the mean droplet size Dj, shown in Eqn. (2). 20 
Values for these four parameters are currently only available for oxygen/hydrogen propellants. 20 A procedure that 
can be used to recalculate these parameters for oxygen/methane propellants is described in Ref. [20]. In the current 
model exercises, these parameters were varied to determine sensitivities of the analysis techniques and determine 
whether other values could best represent the test data. 
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Atomization 

Coefficient 


(i) 


D. 

) 





P /r 


Drop Size 

Coefficient 



2/3 


( 2 ) 


For combustion analyses with ROCCID and CICM, injector mixing plays a key role that is often misunderstood. 
In the distributed energy release calculation procedures described in the JANNAF performance methodology, 25 as 
used in these models, the injection process is treated separately from the combustion process, which means that mass 
and mixture ratio distributions from the injector are defined before combustion is initiated. The part of the thrust 
chamber energy release efficiency associated with mixing loss, 7j mix , is thus defined in the JANNAF performance 
methodology as the product of several independent factors, including the injector inter-element mixing efficiency, 
and the global mass and mixture ratio bias, which can be intentional (such as mass or mixture ratio variations across 
the injector face for combustion stabilization or near the chamber wall for wall compatibility), or unintentional 
(mass and mixture ratio variations across the injector face as a result of maldistributions from manifolding, element 
geometrical tolerances, etc.). In ROCCID, such injector face mass and mixture ratio biases can be included by the 
use of the two-zone (core and barrier), four-streamtube (core, barrier, baffle, and fuel film coolant) inputs. 19 In 
CICM, these zone biases must be calculated with independent runs and mass weighted outside the program, as 
suggested by Refs. [20] and [25]. However, neither program provides a completely model-defined calculation of the 
inter-element mixing efficiency. In the ROCCID model, an zVzfer-element or overall injector mixing efficiency is 
defined by the user. This mfer-element mixing efficiency can be calculated from m^ra-element - or single element - 
Rupe mixing efficiency, 26 or E m , using techniques developed by Ito and Calhoon 25 combined with methods from the 
Liquid Injection Spray Pattern (LISP) calculation 28 or from Nurick. 29 No mixing calculations are included implicitly 
in CICM, which essentially calculates only the vaporization efficiency of the single element streamtube, which for 
single element type patterns can be assumed to be the vaporization efficiency of the thrust chamber, 7j vap , as defined 
in the JANNAF performance methodology. The overall thrust chamber energy release efficiency (the product of 
rj mix and ij vap ) is provided in the ROCCID code, but must be calculated independently when using CICM. 


1 . Single Element Cases - Shear Coax 


Modeling of the single-element shear coaxial element injector data showed how energy release rate (7j c * or rj Isp ) 
profiles, compared to axial and circumferential local heat flux profiles. The intra-element E m for the shear coaxial 
element cases were derived from a cold flow correlation developed by Burick, 30 which is reproduced in Fig. 2. For 
this correlation, E m was measured 5 inches from the injector face with a two-phase deceleration probe in a 
backpressure chamber, with water and air simulants for oxidizer and fuel, respectively. 30 Note that for most of the 
L0 2 /gCH 4 element conditions, the correlation in Fig. 2 is out of range, as indicated. 

A ROCCID calculation for Case 1 from Table I using an input E m of 83.5% and the internal submodel selection 
of DROPMIX droplet size parameters is shown in Fig. 3. For these input conditions, the predicted rj c * matches the 
measured tj c * and all the energy release loss is from mixing and none is from vaporization. Note that for Case 1 
element operating conditions, the Burick correlation for E m is out of range. Extrapolating the Burick correlation to 
Case 1 conditions suggests an E m of about 81%, which would underpredict the measured rj c * by about 0.7%. 

In the single-element experiments, axial and circumferential local heat flux profiles were measured, 9 which may 
be considered a reasonable representation of the energy release rate. Figure 4 shows a 30% increase in droplet size 
from the DROPMIX submodel improves the match of the rj c * profile to the heat flux profile and still allows a match 
of measured and calculated rj c *. Further increases in droplet size result in a vaporization efficiency less than 100% 
and a mismatch between measured and calculated rj c * unless E m is also increased. 
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Figure 2. Correlation of cold flow single-element shear coax injector element Rupe mixing efficiency (£ m ), from 
Burick, 30 extended to L0 2 /gCH 4 . Solid data points - oxidizer post recess of one orifice diameter. Open data 
points - no oxidizer post recess. 



Axial Position, inches 



Figure 3. ROCCID analysis of single element 
recessed shear coax (Case 1). 


Figure 4. ROCCID analysis of single element recessed 
shear coax (Case 1) with increased droplet size. 


Similar calculations with ROCCID for a smaller element without and with a cup recess are shown in Figs. 5 and 
6, respectively. For Case 2 from Table I, the lack of the recessed oxidizer post (“cup”) significantly decreased 
measured ?j c * as well as heat flux near the injector face. The Burick correlation of E m for non-recessed shear coaxial 
element is within the range of Case 2 conditions, and using the predicted E m of 83% resulted in the requirement to 
increase the DROPMIX droplet size by a factor of 2.2, as shown in Fig. 5. The mixing loss was about 2.5% Arj c *, 
and the vaporization loss was about 3.5% A rj c *. For Case 3 from Table I, the same element but with a recessed 
oxidizer post, the Burick correlation of E m for recessed shear coaxial element is also within the range of Case 3 
conditions, and the predicted E m of 84% and the 1.3 factor on DROPMIX used for Case 1 compares favorably with 
measured rj c * and heat flux profiles, as shown in Fig. 6. Use of the anchoring parameters for Case 3 compare very 
well with use of the anchoring parameters for Case 1 , even though Case 3 has a smaller element and about half the 
fuel-to-oxidizer injection velocity ratio. However, the Burick parameter to predict E m is about the same, which 
suggests that the Burick correlation provides a reasonable means to predict mixing for recessed shear coaxial 
elements when in range of the correlation. 

Off-nominal operating conditions for the recessed single-element shear coaxial injector of the single-element 
PSU testing, shown as Cases 9-12 in Table I, were next evaluated. The effect of throttling on tj c * for this and the 
other injectors from the single-element test program was shown in Fig. 1. For the recessed shear coax design 
represented in Cases 9-12, at a MR of 3.0, the performance appears to peak at a P c of 1000 psia. 
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Figure 5. ROCCID analysis of single element non- 
recessed shear coax (Case 2). 


Figure 6. ROCCID analysis of recessed single 
element shear coax (Case 3). 


The results of modeling a low-i^ single-element recessed shear coax case with ROCCID are shown in Fig. 7. 
For this case, because the fuel density is lower, the E m calculated from Burick is in the range of the correlation but 
high - about 89% - and approaching the maximum. However, atomization quality has worsened, so the 
improvement in mixing efficiency is offset by a decrease in vaporization efficiency. Still, a combination of an E m of 
89% and a 30% increase to the droplet size calculated by the DROPMIX program significantly overpredicts the 
performance. The performance mismatch can be corrected by using a 2.4 factor increase in droplet size with the E m 
of 89%, as shown in Fig. 7(a), or an E m of 79.5% and the 1.3 factor increase to the droplet size, as shown in Fig. 
7(b). The former appears a better match to the heat flux data profile. 



0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

Axial Position, inches 


(a) E m = 89% and 2.4x increase in dropsize 


o 


2 3 4 5 6 7 8 9 10 11 12 13 14 15 

Axial Position, inches 



(b) E m =19.5% and 1.3x increase in dropsize 


Figure 7. ROCCID analysis of single-element recessed shear coax injector at 7^-300 psia (Case 9). 


The results of modeling a high-7^. single-element recessed shear coax case (Case 12) with ROCCID are shown in 
Fig. 8. For this case, because the fuel density is much higher, an E m from Burick is considerably out of range of the 
correlation; an estimate for E m is 76%. Thus, mixing efficiency is considerably decreased from the low-7^. case, 
although the %* loss has been partially offset by an improvement in atomization quality and increase in vaporization 
efficiency. Still, with use of an E m of 76% the performance is underpredicted by 1.5% even with 100% vaporization 
efficiency. Use of an E m of 81% allows the calculation to match the measured data. 

Thus, modeling off-nominal conditions of the single-element recessed shear coaxial injector with ROCCID 
suggests that using the Burick correlation for E m - either reading directly from the plot or extrapolating - may be 
erroneous, especially at low-7^ conditions. Using a 30% increase to the droplet size calculated by the DROPMIX 
program, E m is required to be about the same value — 82% - for all cases. A re-evaluation of the mixing efficiency 
from Burick for elements with dense fuels is recommended. 

CICM is currently capable of analyzing only shear coaxial injection elements with liquid propellant in the center 
jet and gaseous propellant in the annulus, and has only been validated with liquid oxygen and gaseous hydrogen 
propellants. 20 These restrictions makes the model less appealing for developing capability to analyze L0 2 /LCH 4 , but 
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Figure 8. ROCCID analysis of single-element recessed shear coax injector at 7^-1200 psia (Case 12). 

it was included in this analysis development program since CICM is the industry standard for liquid/gas shear 
coaxial element calculations in JANNAF, 25 and its calculation procedure for atomization may be used in future 
models. Note that as discussed earlier, CICM does not calculate or include an intra-element mixing non-uniformity, 
such as shown in Fig. 2, so that the mixing loss must be included either through averaging of several calculations at 
different mixture ratios, or an ad hoc addition of a mixing efficiency to the CICM-calculated vaporization 
efficiency. Such mixing losses were not included in the results presented here. 

The use of CICM for oxygen and methane calculations requires new propellant and combustion property data, 
and new parameters for the semi-empirical atomization correlations. The thermodynamic property data required for 
CH 4 propellant and L0 2 /CH 4 combustion have been developed during this program. Atomization parameters are 
currently only validated in available documentation for oxygen/hydrogen propellants with shear coaxial element 
geometries similar to those in the J-2 injector. 20 Three methods to recalculate these parameters for other propellants 
or radically different shear coaxial element geometries, which all use the cup pressure drop in the calculations, are 
described in the CICM User’s Manual. 20 Unfortunately, the methods were found to be inapplicable to most of the 
single-element L0 2 /gCH 4 experimental test cases since the cup pressure drops were found to be practically zero, 
which is not surprising given that reaction rates for L0 2 /gCH 4 propellants are considerably slower than for L0 2 /gH 2 
propellants. The best opportunity to determine cup pressure drops was with the elements in Cases 2 and 3, using 
both recessed and flush configurations. Unfortunately, even for this case, the cup pressure drop could not be 
determined for several cases of P c and MR. Thus, the recommended methods in Ref. [20] to determine new 
atomization parameters could not be used. 

Consequently, to determine correlating parameters for the L0 2 /gCH 4 propellants, all of the parameters were 
varied individually, starting with the values in the CICM user manual for L0 2 /gH 2 propellants. Two parameters 
were found that showed the strongest influence and provided the most correct trends - chamber droplet formation 
size parameter BSPRC, and chamber liquid jet stripping rate CSPRC. These two parameters were subsequently 
varied together, in an attempt to determine the best match of the heat flux profile and the r/c* data from all the 
appropriate single-element test cases. All of the cases were then re-examined with the best correlating parameters. 

The initial calculation using manual-suggested inputs for B A and C A in the cup and the chamber from L0 2 /gH 2 
propellants is shown in Fig. 9. Comparing the heat flux profile and the profile calculated with these parameters 
shows the energy release near the injector (the first 2 or 3 inches of the combustion chamber) may be 
underpredicted, while the final calculated %* significantly underpredicts the measured ri c * even though the result is 
from vaporization losses only and does not include mixing losses as found in a parameter such as E m . Use of the B A 
and C A in the cup and the chamber from the manual are obviously not valid with oxygen and methane propellants. 

Selection of better values for atomization parameters was subsequently an exercise trading improving the 
atomization rates and hence increasing the vaporization rates of the cases with under-predicted performance, and not 
significantly altering the already higher performing cases. A prediction with atomization parameters revised for 
L0 2 /gCH 4 is also shown in Fig. 9, which indicates a more rapid rise in energy release rate near the injector. The 
chamber liquid jet stripping parameter had the largest effect on the initial vaporization rate near the head end of the 
chamber, while the chamber droplet formation size parameter was used to reduce the vaporization rate in the far 
field. However, while there was a large change in the profiles for the element in Case 1, there was only a minor 
change to the profile for the smaller geometry element in Case 3. The revised atomization parameters provided 
similar modeling comparisons for most of the cases analyzed. 
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Figure 9. Comparison of CICM atomization parameters 
to match axial rj c * profiles to axial heat flux data. 
Case 1 at TVT000 psia. 


Figure 10. CICM L0 2 /gCH 4 atomization 
parameters underpredict 7j c * for Case 9 at 
,P C ~300 psia. 


However, using the revised atomization parameters for low TV-300 psia (Case 9 in Table I), CICM significantly 
underpredicted rj c *, as shown in Fig. 10. One possible reason for this error is that oxygen is subcritical at this P c 
compared to other tests examined. Correlation parameters from conditions with much lower (subcritical) chamber 
pressures are likely to be different than those developed for high chamber pressures. Second, this P c is throttled to 
30% of the nominal design point of this injector element, providing significantly off-nominal operating conditions. 

2. Single Element Cases - Swirl Coax 

For the swirl coaxial single-element injector, Case 4 from Table I, several inputs had to be modified for the 
measured rj c * and the axial heat flux profiles to match. Results from unmodified calculations are shown in Fig. 11. 
The E m for the swirl coaxial elements were estimated from data in Hulka et.al., 31 for which E m was measured 5 to 6 
inches from the injector face with a patternator at atmospheric pressure, with water and sucrose simulants for 
oxidizer and fuel, respectively. Swirl coaxial element geometries were varied at constant oxidizer-to-fuel 
momentum ratios of about 3, and found the oxidizer free swirl angle to be most influential on E m ? x Figure 1 1 shows 
the r/c* was overpredicted by 1.2%, while the start of combustion is delayed too far downstream due to a significant 
overprediction of primary atomization lengths. To improve the match to measured parameters, the program- 
calculated atomization length had to be reduced by a factor of 10, the droplet size from the AEROJET submodel 
increased by a factor of 2.7, and the E m estimated from Ref. [31] reduced from 94% to 88%. The results of these 
modifications are shown in Fig. 12. A large change in E m is not surprising given the differences in oxidizer-to-fuel 
momentum ratio for Case 4 (about 0.3) from Ref. [31] (about 3). The significant alterations from calculated 
atomization parameters (primary atomization lengths, droplet diameters used in vaporization calculations) suggest 
that the atomization model does not adequately represent the mechanisms of the swirled oxidizer atomization. 




Figure 11. ROCCID analysis of single-element swirl 
coax (Case 4), without modifications. 


Figure 12. ROCCID analysis of single-element swirl 
coax (Case 4), with modifications. 
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3. Multielement Cases 


Calculations with ROCCID for a multielement liquid/gas shear coaxial injector 7 are shown in Fig. 13, using the 
same 1.3 factor on DROPMIX droplet size from the previous calculations with recessed shear coaxial elements, 
shown in Figs. 4 and 6. However, these conditions are far outside the range of the Burick correlation of E m for 
recessed shear coaxial elements shown in Fig. 2. An extrapolation of the Burick correlation to these conditions 
would result in a completely unrealistic E m , so a prediction of E m could not be made. An E m of 83.6% was used to 
match measured rj c *. Clearly, if shear coaxial injection elements will be used for future injectors using liquid 
oxygen and gaseous methane propellants, the Burick correlation for E m should be extended to include typical ranges 
of values for those propellants. While there were no heat flux data for this element, CICM analysis with the 
L0 2 /CH 4 correlating parameters provides a reasonable match with the ROCCID results, as shown in Fig. 14. 



Axial Position, inches 



Distance from I njector Face, inches 


Figure 13. ROCCID analysis of multi-element shear 
coax injector. 


Figure 14. Comparison of CICM and ROCCID tj c * 
profiles for multi-element shear coaxial injector. 


B. Combustion Stability 

With one exception, combustion instability with L0 2 and gCH 4 propellants has been relatively absent. In one 
82-element shear coaxial element injector in a 5.66” combustion chamber diameter at booster-level P c ~\ 500 psia to 
2000 psia, injection-coupled combustion instability occurred as a result of coupling between the “organ-pipe” 
resonance of the oxidizer element tube and the first tangential (IT) mode of the combustion chamber. 7 The 
comparison of organ pipe responses for injectors with different length injector oxidizer tubes and the combustion 
chamber response is shown in Fig. 15. Further details and discussion of the injection-coupled analyses of this 
instability are provided in Refs. [7] and [32]. 



(a) Injector configuration that exhibited spontaneous 
IT mode combustion instability 



Fitqumtf r Mi 

(b) Injector configuration that was stable 


Figure 15. Comparison of calculated combustion chamber and oxidizer injector injection responses for 
two shear coaxial element injectors with different oxidizer post lengths. 7 
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C. Regenerative Chamber Heat Transfer and Cooling 

In engine cycles with regeneratively cooled combustion chambers or nozzles, the methane delivered to the 
injector is in gaseous form because of heat addition as a coolant. Development of such chambers will require further 
research and development into how subcritical methane changes phase in the coolant circuits, and especially how 
that process couples to other engine dynamics and becomes a system oscillation. Additional development of the 
methane combustion chamber cooling design features is also needed. NASA is currently conducting work in these 
areas, with plans to report in future papers. 


III. Liquid Oxygen and Liquid Methane 


The use of pressure-fed engine systems, such as expected for in-space applications, usually means that 
combustion chamber or nozzle regenerative cooling is not likely or not feasible, so the methane will be delivered to 
the injector in liquid form. With both propellants now in liquid form, the impinging element geometry, from 
storable propellant or oxygen/RP-1 history, can be employed. Reasonable combustion efficiencies have been 
demonstrated along with stable combustion, although modeling has suggested that higher efficiencies can be 
achieved. However, because the fuel density is not as high as MMH or RP- 1 , the coaxial element geometry can also 
be considered. In-house NASA MSFC testing with swirl coaxial element injectors in 5.66” diameter combustor 
hardware has demonstrated high efficiencies in relatively long combustion chambers (exceeding 20 inches in 
length), although combustion instability, at both low and high frequencies, has occurred with several 
configurations. 15,17,18 The coaxial element used in these studies swirled the L0 2 in the center post and injected the 
LCH 4 axially through the outer annulus, as shown in Fig. 16. The oxidizer metering orifice, the fuel sleeve, and the 
face nut were modular and replaceable allowing for variations of oxidizer and fuel metering area to be tested. The 
tangential entry orifices for the oxidizer swirl were slots formed by electro-discharge machining. 


Fuel sleeve Fuel facenut 



Figure 16. Typical swirl coaxial injector element details. 


15,17,18 


A. Performance 

Performance calculations of tests using L0 2 /LCH 4 propellants with a NASA MSFC in-house multi-element like- 
on-like (LOL) doublet impinging element injector 15 were conducted with ROCCID. E m for this element type can be 
predicted from several sources, including Refs. [28], [29], or [33]. Using the like-impinging doublet correlation 
from Ref. [29] with no edge-on-edge offset and the design cant angle, an E m of 77% is predicted. Using this E m and 
the DROPMIX atomization parameters, the initial calculation is shown in Fig. 17. The predicted rj c * of 95.4% in the 
initial model without modifications provided an excellent match to the measured rj c * of 95.6%. 
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Figure 17. ROCCID analysis of multi-element paired like-on-like injector with L0 2 /LCH 4 propellants. 
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The good match of ROCCID predictions to measured performance with this injector concept suggests several 
avenues to increase performance. First, several aspects of the paired doublet geometry can be altered to increase the 
E m , e.g., an increase in the cant angle. Within the range of the correlation from Ref. [29], E m can be increased from 
77% to 84% and rj c * can improve up to 2%. References [28] and [29] also suggest impingement angles and 
momentum ratios are not optimum for performance. Second, the orifice dimensions could be reduced and the 
number of elements increased, according to LOL element combustion stability guidelines such as from the Hewitt 
Stability Correlation, 34,35 and still meet combustion stability margin requirements, especially if baffles and/or head- 
end acoustic resonator cavities are included. While the tested element packaging had an element pattern density 
similar to several development and flight L0 2 /RP-1 LOL element injector patterns, higher performance (-99% 77 c *) 
has been achieved in drilled-orifice element patterns with 40% more elements-per-area, such as the Titan I second 
stage engine. Higher element density patterns can also be fabricated with other means than conventionally drilled 
plates, such as platelets or micro-orifices. The use of an improved E m and more injection elements suggests the 
performance of this injector can be significantly improved. On the other hand, a wall compatible element or some 
level of barrier cooling would subsequently be required which would reduce performance. 

Calculations of a test from a NASA MSFC in-house 28-element swirl coaxial element injector 15,17 are shown in 
Fig. 18. To be able to perform these analyses, the ROCCID code was modified in 2008 to include the capability to 
analyze a swirl coaxial injector element with liquid flow from both circuits. 19 In Fig. 18, the uni-element E m - 94% 
was estimated from Hulka et.al., 31 a cold flow mixing study conducted with liquid/liquid simulants. However, the 
new atomization submodel provided unreasonably large atomization lengths and droplet sizes that resulted in greatly 
underpredicted %*. For the predictions to match the test data, considerable reductions were required in atomization 
length (factor of 0.2) and droplet size from the AEROJET model (factor of 0.17), similar to alterations shown earlier 
from the L0 2 /gCH 4 cases. A calculation using the same atomization factors, at similar P c and MR but with a longer 
combustion chamber, is shown in Fig. 19. Since the measureable increase in rj c * from the extra 4” of combustion 
chamber length was predicted reasonably well, the injector is shown to be vaporization limited, and the droplet size 
parameters used in the calculations are probably representative. Note that the L0 2 vaporization efficiency is lower 
than the LCH 4 vaporization efficiency, with both contributing to performance loss in the L -20” chamber. At 
L -24”, the LCH 4 vaporization is practically complete, while the L0 2 vaporization is still limiting performance. 




Axial Position, inches 


Axial Position, inches 


Figure 18. ROCCID analysis of 28-element swirl 
coax injector with L0 2 /LCH 4 propellants. 


Figure 19. ROCCID analysis of 28-element swirl 
coax injector with L0 2 /LCH 4 propellants. 


Tests with a 40-element injector with swirl coaxial elements and liquid/liquid propellants where the fuel injector 
flow paths were redesigned showed increased performance (and apparent reduced chug oscillation amplitudes, 
although high-frequency pressures were not measured in this phase of the test program). 18 Modeling of these cases 
is slightly complicated by the diversion of some fuel from the injector element fuel flow to fuel film coolant (FFC) 
at the chamber wall. Consequently, in the ROCCID input file the injector pattern was split into a core zone of the 
inner rows and a barrier zone of the outer row into which the FFC was mixed. To provide reasonable performance 
calculations, ROCCID needs to mix the FFC into a barrier zone rather than leave it unmixed in an outside 
streamtube. 

The effect on performance by altering the fuel injector coaxial annulus (or “gap”) width can be examined with 
test data from the 40-element swirl coax element injector. 18 rj c * improved from -90-91% to -95-96% with a 
reduction in fuel gap width. 18 Modeling with ROCCID for this geometrical variation was attempted but the 
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liquid/liquid swirl coaxial element atomization model failed to produce rational values of either fuel or oxidizer 
droplet sizes for the larger fuel gap widths, and reasonable values for the smaller fuel gap widths. User-defined 
droplet sizes were used to run the model, but little mechanistic understanding could otherwise be derived about the 
droplet size. The observed improvement in from reducing the fuel gap widths can come from several sources - 
including improvement in vaporization efficiency due to better fuel atomization, improvement in vaporization 
efficiency due to better oxidizer atomization, and improvement in mixing efficiency due to better unielement 
mixing. While all possibly contribute to the performance improvement, the most likely is improvement in 
vaporization efficiency due to improved fuel and oxidizer atomization. An example of a ROCCID output from a test 
with a larger fuel gap width, and which includes heat flux data from a calorimeter chamber, is shown in Fig. 20. 
This calculation uses very high mixing levels as suggested from Ref. [31], but relatively large mass median droplet 
sizes. It is not obvious whether that is a good assumption. However, comparison of the axial rj c * profile to the heat 
flux data suggests the selected model parameters are reasonable. 



0 5 10 15 20 25 


Axial Position, inches 


Figure 20. ROCCID analysis of larger fuel gap width 40-element swirl coax injector with 
L0 2 /LCH 4 propellants. 


A calculation of a test of a 40-element liquid/liquid swirl coaxial element injector with a smaller gap width in an 
L -21” chamber is shown in Fig. 21. The predicted value of uni-element E m ~ 95% from Ref. [31] is used, along 
with considerable reductions in atomization length (factor of 0.2) and droplet size from the AEROJET models 
(factor of 0.155). A calculation with the longer combustion chamber ( L -25”) is shown in Fig, 22. Similar to Figs. 
18 and 19, rj c * continues to improve with the increase in chamber length, suggesting the 40-element injector with the 
small fuel gap width is still vaporization limited. For these 40-element injectors, another possible performance 
improvement with these injector tests may be increased entrainment of the FFC into the core (i.e., improvement of 
inter- streamtube mixing). However, such entrainment is typically slow, and very unlikely to change sufficiently 
between 21” and 25” downstream of the injector face to increase t] c * by the margin observed. 



Axial Position, inches 

Figure 21. ROCCID analysis of 40-element swirl 
coax injector with smaller fuel gap width and 
L0 2 /LCH 4 propellants. L ~21”. 



Axial Position, inches 

Figure 22. ROCCID analysis of 40-element swirl 
coax injector with smaller fuel gap width and 
LO 2 /LCH 4 propellants. L ~25”. 
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B. Combustion Stability 

1. Low Frequency 


Nearly every test with a NASA MSFC in-house 28-element swirl coax element injector using L0 2 /LCH 4 
propellants experienced some level of low-frequency combustion instability, or “chug.” 15,17 However, while the 
chug frequency was dominant in these tests, the waveforms of the chug oscillations often included significant high- 
frequency content, especially near the maximum pressure of the chug oscillation cycle. Typical waveforms from a 
lower P c test are shown in Fig. 23. Usually, a chug waveform is reasonably sinusoidal and not nearly so 
contaminated with high frequency signals. Counting higher frequency waves in Fig. 23b suggests the higher 
frequency is about 4000 Hz. A frequency spectrum of this signal is shown in Fig. 24, where 4000 Hz is also 
identified by the acoustic-like rounded signature. The spectrogram is calculated from 32768 samples over the time 
window as shown, with no filtering applied. 



(a) View showing low-frequency content 



Time from Start, seconds 


(b) Close-up showing high-frequency content 


Figure 23. Chug waveform during Test P25 14-042 with 28-element injector. 
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Figure 24. Spectrogram showing chug during Test P25 14-042 with 28-element injector. 

Throttling up (increasing the P c ) and lowering the MR significantly reduced the chug amplitude. High-frequency 
content was still present, but the frequency characteristics at this higher P c were different than from tests at the lower 
P c , as will be described in a later section. In a combustion chamber approximately with half the length, significantly 
increased chug amplitude was observed, while in a longer combustion chamber chug amplitude was reduced. Both 
changes were predictable from changes in L *. 

Similar chug characteristics except with higher overall amplitudes were encountered during tests with an initial 
40-element swirl coax element injector. 17 Following this test, several broken oxidizer swirl metering orifices were 
found loose in the oxidizer manifold, and several other metering orifices were found with bulged or shattered slots. 17 


: Te 
Tir 

st P251 4-042 

— Pc - P0103 

neslice: 5.43 to 5 75 sec 



^ chug at 208 Hz 

















-ill 



1 

i 




14 

American Institute of Aeronautics and Astronautics 


An injector element redesign which reduced the oxidizer slot metering area and revised the fuel sleeve geometry 
also experienced chug, and, later in the same test, spontaneous IT combustion instability. 17 The higher oxidizer 
pressure drop reduced the chug amplitude, while the character of the high-frequency content in the spectrum 
appeared unchanged. Following this test, several broken oxidizer swirl metering orifices were again found loose in 
the oxidizer manifold, with another metering orifice also found with bulged slots. 17 

Initial tests with a further redesigned 40-element injector pattern using the reduced oxidizer metering 
configuration that had previously been unstable at the IT mode, plus a change of the fuel sleeve design, showed 
similar or worse chug characteristics during initial tests. 18 Unfortunately, high-frequency pressure was not measured 
during these initial tests, but only during a limited test series with a calorimeter chamber. Later variations of the fuel 
annulus design with reduction of the fuel gap widths in later tests appeared to improve chug stability characteristics. 
However, no high-frequency pressure was measured in these tests, and static pressure measurements recorded on 
high-speed data systems appeared unreliable as several tests had inconsistent results - i.e., considerable chug 
amplitude on one test and then little on a later test with similar test conditions. These data are still under review. 

Chug test data were compiled and compared to attempt to determine the experimental neutral chug stability 
curve, the measured boundary between stable and unstable operation. Compiled data for the 28-element injector are 
shown in Fig. 25. Data are either marginal (M) or unstable (U/S) except for one test, and would be at or below the 
neutral stability curve. Unfortunately, a second test series using the same injector also had a 20% reduction in fuel 
injector admittance (flowrate/square root of pressure drop), which compromises the analyses. The oxidizer injector 
admittance was similar between the test programs, so the chamber pressure during chug was not a suspected error. 
No explanation for the change in the fuel injector admittance has been found. 17 
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Figure 25. Compiled chug data and model prediction for the 28-element injector. Data = points, 
model prediction = lines. 


Modeling of the chug was conducted with the techniques described by Wenzel & Szuch. 36 This analysis 
procedure defines the neutral chug stability curve from the solution of a linearized characteristic equation describing 
instability in a lumped-parameter impedance combustion chamber. Each propellant has its own specific time lag, a 
delay time which is a function of injection thermodynamics, chamber pressure, and mass flow rates. A mixing time 
lag common to both propellants is also included. Errors in the characteristic equation shown as equation (31) in Ref. 
[36] were corrected in the current analyses. The model was applied to match the characteristics of the neutral 
stability curve defined from evaluation of the data in Fig. 25, including the observed frequencies of oscillation at the 
respective chamber pressures. The results of the modeling are shown in Fig. 25. A solid black line represents the 
neutral stability curve defining the insipient chamber pressure at chug, and a red line represents the frequencies at 
the neutral stability curve. 

To force the model to generate the predictions shown in Fig. 25, only a portion of the total measured fuel injector 
pressure drop was allowed, which means some of the pressure drop was not participating in resisting chug. The fuel 
annulus sleeve designs for the data in Fig. 25 were unusual for coaxial injector design in that most of the fuel 

15 

American Institute of Aeronautics and Astronautics 


injector pressure drop was located in the fuel sleeve holes rather than in the annulus or exit dynamic head. 15,17 
Unlike for liquid oxidizer/gaseous fuel injectors, increased fuel pressure drop for liquid oxidizer/liquid fuel injectors 
is stabilizing. However, to stabilize the system shown in Fig. 25 for expected conditions (about P c = 250 psia) will 
require increasing AP/P C for both injector circuits. 

A revised fuel annulus sleeve design that applied the majority of fuel injector pressure drop to the annulus or exit 
dynamic head rather than the fuel sleeve holes was used in the later tests with the 40-element injector, along with a 
higher-pressure drop oxidizer circuit. 18 Unfortunately, for these tests high-frequency pressure was not measured, so 
reliable dynamic data were not available. Preliminary evaluations with static pressure recorded at high sample rate 
suggest that increased LCH 4 injector pressure drops at the injector element exit and increased L0 2 injector pressure 
drop improved chug stability margin for these tests. Final analysis results of these data are not yet available. 

2. High Frequency 

No test with the 28-element swirl coax element injector exhibited sustained combustion instability at a 
longitudinal or transverse acoustic mode of the combustion chamber. However, high-frequency content, especially 
between 2500-4000 Hz, was always present in the chamber pressure, and contaminated the low frequency signal 
during chug oscillations, as shown for example in Fig. 23b. If the chug frequency would be removed by high-pass 
filtering, the remaining higher frequency waveforms in Fig. 23b would resemble resurge instability. Whether what 
is seen in Fig. 23b is resurge, especially when compared to other 28-element injector tests, is still under 
investigation. More detailed stability data and analyses from these tests are available in Ref. [37]. 

Two different forms of high-frequency content were observed during tests with the 28-element injector. The 
first, at lower P c s (270-300 psia), had a single rounded response with a frequency at the peak amplitude around 4000 
Hz, and other minor content except for 1100 Hz. The second, at higher P c s (315-375 psia), had several distinct 
peaks centered around two frequency ranges - one at 2400 Hz to 3600 Hz, the other at 4800 Hz to 5300 Hz. Planar 
frequency intensity plots for representative tests of each form are shown in Figs. 26 and 27, and spectrograms for 
these tests are shown in Figs. 28 and 29. 

For the test shown in Fig. 28, an equilibrium IT mode acoustic frequency is calculated to be about 5015 Hz, and 
a first longitudinal (1L) mode acoustic frequency is calculated to be about 1240 Hz. Figure 28 shows two obvious 
high-frequency signatures at about 1100 Hz and 4000 Hz, and possibly several others including 2800 Hz, 3500 Hz, 
and 5500 Hz. Longitudinal and transverse acoustic modes that are not unstable are often observed at lower 
frequencies than calculated by the equilibrium combustion sound speed, so the observed frequency of 4000 Hz (or 
an effective sound speed -80% of equilibrium) would not be unrepresentative for a IT acoustic mode. Likewise, at 
an effective sound speed -90% of equilibrium, the 1 100 Hz signal might be the 1L mode. Unfortunately, there was 
insufficient high-frequency pressure instrumentation in the combustion chamber to verify that what was observed at 
4000 Hz was the IT mode. 

The RMS pressure amplitudes for the signal at 4000 Hz shown in Fig. 28 represent fairly small amplitude waves, 
typical for an unexcited acoustic mode. However, the amplitudes of the 4000 Hz oscillation at the maximum 
pressure of the chug wave as shown in Fig. 23(b) are much larger, - 20 psi peak-to-peak (-7 psi RMS), or almost as 
large as the amplitude of the chug oscillation itself. The spectrograms in Figs. 24 or 28 thus are not good portrayals 
of the real character of the high-frequency oscillations shown in Fig. 23. 

The spectrogram for the higher P c test shown in Fig. 29 has a different character than for the low P c test shown 
in Fig. 28. Equilibrium acoustic frequencies were nearly identical with the lower P c tests, despite considerable 
differences in P c and sometimes MR. Figure 29 identifies four or five peaks of acoustic-like behavior at the higher 
P c - 5100-5200 Hz, 3400-3500 Hz, 2500-2800 Hz, and 1000 Hz. The chug is still identifiable between 165 and 230 
Hz. Closer examination of raw waveforms of several of these higher-pressure 28-element injector tests reveals that 
while there was no sustained high-frequency instability, several tests displayed bursts of short-duration spontaneous 
instability. Results of an instructive test, at high P c and low MR where peak-to-peak chug amplitudes had been 
reduced to less than 1% of mean P C9 is shown in Fig. 30. In this test, numerous periods of extended spontaneous 
high-frequency instability occurred, lasting 15 to 40 msec, several indicated by the arrows in Fig. 30. A close-up of 
a waveform during one of those bursts is shown in Fig. 31, clearly indicating multiple frequencies. Figure 32 
compares spectrograms of a narrow time slice of one of the bursts to a longer average time. The spectrogram over 
the longer time slice is similar to Fig. 29, showing peaks of high-frequency behavior at 5100 Hz, 3400 Hz, 2500 Hz, 
and 990 Hz. However, the spectrogram over the narrow time slice in Fig. 32 shows that two of those signals - at 
2540 and 5080 Hz - have been excited, along with another signal at 7620 Hz. The 5080 Hz signal in Fig. 32 could 
be the IT mode, nearly exactly at the equilibrium value, which may have been excited by the signal at 2540 Hz. 
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Figure 26. Planar spectrogram of high frequency pressure from low-,P c test with 28-element 

injector. Test duration increases from top to bottom. Color scale at right is decreasing RMS 
pressure from top to bottom. The plot at left shows static chamber pressure. 
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Figure 27. Planar spectrogram of high frequency pressure from high-,P c test with 28-element 
injector. Test duration increases from top to bottom. Color scale at right is decreasing RMS 
pressure from top to bottom. The plot at left shows static chamber pressure. 
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Figure 28. Spectrogram of low P c test with 28- 
element injector. 
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Figure 29. Spectrogram of high P c test with 28- 
element injector. 
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Figure 30. Pressure time trace during Test P2514- 
045 with 28-element injector. 
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Figure 31. Waveform of high-frequency activity 
during Test P25 14-045. 
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Figure 32. Spectrograms from Test P25 14-045 during high-frequency activity. 


Since the bursts of high-frequency activity shown in Figs. 30 to 32 were not sustained, and were not the result of 
dynamic stability rating perturbations, they present some challenges for categorization when consulting combustion 
stability guidelines such as Ref. [38]. The oscillation at 5.08 seconds in Fig. 30 is obviously discrete and well above 
the amplitudes of the noise floor or the averages of the transverse acoustic modes. The peak-to-peak amplitude of 
the oscillation at 5.08 seconds is 3% to 5% of the mean chamber pressure, relatively small but not typically 
acceptable for transverse mode stability for flight vehicles. Had the high-frequency event at 5.08 seconds been the 
result of a dynamic stability rating perturbation, the damp time of 37 msec would have exceeded the allowable damp 
time according to Ref. [38], calculated to be about 18 msec for the frequency at 5080 Hz if that is the IT mode. 

Understanding the high-frequency combustion stability characteristics of these tests requires identifying the 
sources of the signals shown on Figs. 26-32. Because of insufficient high-frequency pressure instrumentation during 
these tests, it is not readily apparent what these signals represent. Sources of such signals can include sensor tube 
resonances, dynamics of injector element features such as the oxidizer tube or the fuel annulus, or dynamics of 
combustion processes such as atomization, vaporization, mixing, and reaction. Two can be quickly excluded: 1) the 
sensor tube resonance of the high-frequency pressure transducer, according to sensor port configurations, is about 50 
kHz, and 2) the quarter- wave resonance of the fuel annulus of the injector element is about 1 1 kHz. 

Axial (or “organ pipe”) resonances of the flow in the oxidizer tube of a shear coaxial injector element have been 
recognized for some time as sources of combustion instability in coaxial element injectors, 7 ’ 32,39,40 and calculating 
such resonances in a shear coax element oxidizer tube is fairly straightforward. However, the dynamics of a hollow- 
core swirl flow in a tube are more complicated. 41 As first order estimates, simple calculations of quarter-wave 
(closed/open acoustic) resonances are usually close to the relevant frequencies, since the upstream end of the 
injector element is usually metered and approximates a closed acoustic boundary. For the 28-element injector swirl 
coax element, a quarter- wave resonance of the annular oxidizer film attached to the tube wall could be about 1500 
Hz, while a quarter-wave resonance of the gaseous fluid in the center depends on the composition of the gas and can 
range from 400 Hz with cold gaseous oxygen to 2600 Hz with equilibrium combustion products. These calculations 
suggest the 2500-2800 Hz signals in Figs. 26-32 may be resonances of the gaseous core in the oxidizer tube. 
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Another source for high-frequency combustion instability has been suggested by Hewitt. 34,35 The Hewitt 
Stability Correlation is based on a parameter d/V - injection orifice dimension d divided by injection velocity V of 
the least volatile propellant - that describes observed combustion instability of like-on-like doublets with storable 
propellants. 34 The correlation has been extended to other propellant combinations and element types 42,43 including 
swirl coaxial elements with liquid oxygen and dense gaseous hydrogen. 44 For the latter, with the assumption that 
this element responds like an unlike doublet (due to direct impingement of unlike propellant streams), the injection 
orifice dimension d has been expressed as the film thickness of the hollow-cone swirl flow as calculated by a 
method such as Doumas and Laster. 45 

While in most applications the fuel is the least volatile propellant, for L0 2 /LCH 4 propellants the L0 2 is the lesser 
volatile, although the volatilities are much more similar compared to other propellant combinations such as 
oxygen/hydrogen or oxygen/kerosene. Consequently, for this paper both propellants will be analyzed. For the 
LCH 4 , the injection orifice dimension d is suggested to be either the fuel annulus “gap width” or the hydraulic 
diameter of the fuel annulus, which is the twice the gap width. 

The Hewitt Stability Correlation has been presented two ways, with d/V as a function of chamber diameter in 
which instability was observed, or with d/V as a function of the maximum sustainable frequency of instability which 
can be supported by that d/V , which is a function of element type but not propellant combination. 34,35 Using the 
latter method for both L0 2 swirl film and LCH 4 annular flows predicts a range of maximum sustainable frequencies 
of instability squarely in the range of frequencies observed in the test data. Maximum sustainable frequencies of 
instability from d/V range from 2700-3700 Hz for the L0 2 swirl film (assuming film thickness for d ), and from 
2000-3000 Hz for the LCH 4 annular flow (assuming gap width for d). Similar frequencies from periodic spray 
oscillations have been measured in cold flow tests of hollow-cone swirl flows. 46 

However, if the combustion system is responding to d/V - generated frequencies, the frequencies should show a 
functional relationship to V. Measured frequencies of signals obviously above the noise floor in five frequency 
ranges - around 1000 Hz, 2500 Hz, 2800 Hz, 3500 Hz, and 5000-5500 Hz - were essentially constant versus 
calculated oxidizer and fuel injection velocities. These data suggest that these particular frequencies are likely to be 
resonances of injection element features or acoustic modes of the combustion chamber. Comparing to predictions 
using equilibrium combustion product sound speed, the signals at 1000 Hz are likely the 1L mode and the signals at 
5000-5500 Hz are likely the IT mode. The 2500 Hz, 2800 Hz, and 3500 Hz are thus likely to be resonances of the 
oxidizer post. Certainly a calculation of a 2500 Hz resonance frequency of the swirl coaxial element from a gaseous 
core with combustion product sound speed is encouraging. 

While these discussions provide preliminary rationale for most of the observed frequencies, there is still no clear 
rationale for one of the major differences between Fig. 26 and Fig. 27, which was the disappearance of the 4000-Hz 
signal at lower chamber pressures when throttling up to higher chamber pressures. A rationale for that may be 
provided by examination of 40-element swirl coax element injector test data. 

Three tests with two different configurations of 40-element injectors had high-frequency instabilities near 4500 
Hz apparently at the first tangential (IT) acoustic mode of the combustion chamber. With only one high-frequency 
pressure sensor in the combustion chamber, positive identification of the unstable mode was not possible. However, 
the measured frequency was found to be about 90% of that calculated with equilibrium sound speed for the IT mode 
in all three cases. Also, in one test which occurred in a calorimeter chamber, there was a dramatic increase (about a 
factor of ten near the injector) in measured circumferential heat flux in the combustion chamber. Had the high- 
frequency oscillation been only due to excitation of element resonances, it is unlikely the heat flux would have 
increased so dramatically. These factors suggest that the 4500-Hz oscillation was an unstable IT mode. 
Unfortunately, most of the tests with variations of 40-element injector features did not contain high frequency 
pressure sensor measurements, either in the chamber or the injector manifolds, or any accelerometer data, so it was 
not possible to characterize the high-frequency stability of these tests or their hardware changes. 

The first two apparent IT instabilities with the 40-element injector occurred with the first two tests of a new 
oxidizer element with a smaller metering orifice than used in initial testing. A planar frequency intensity plot of one 
test is shown in Fig. 33. The onset of instability is clearly seen at 4550 Hz, along with the appearance of a harmonic 
at 9100 Hz; higher frequency harmonics can be seen in the spectrum up to a test facility low-pass filter cut-off at 20 
kHz. Waveforms immediately before and during the transition to the instability are shown in Fig. 34, illustrating the 
transition from chug to high-frequency instability. Figure 34(a) shows the spontaneous excitation of a 3910 Hz 
oscillation from one chug cycle, lasting about 8 msec. Note in Fig. 34(b) that the instability was initiated at the 
maximum pressure of the chug oscillation cycle, where similar high-frequency dynamics have been observed from 
nearly all tests. This result resembles the unstable outcome from a perturbation of a dynamic stability rating test. A 
waveform of the high-frequency instability is shown in Fig. 35. Spectrograms immediately before and after the 
instability are shown in Figs. 36 and 37. 
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Figure 33. Planar spectrogram from Test P7061-014 showing high-frequency instability. Test 
duration increases from top to bottom. Color scale at right is decreasing RMS pressure from 
top to bottom. The plot at left shows static chamber pressure. 
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Figure 34. Waveforms from Test P7061-014 of chug and transition to high-frequency instability. 
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Figure 35. Waveforms from Test P7061-014 of high-frequency instability 
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Figure 36. Spectrogram showing low-level chug 
prior to high-frequency instability during Test 
P7061-014 with 40-element injector. 


Figure 37. Spectrogram showing high-frequency 

instability during Test P7061-014 with 40-element 
injector. 


For the high-frequency instability shown in Figs. 33, 35, and 37, peak-to-peak amplitude is observed to be about 
1 0% of the mean chamber pressure, which is not especially large for transverse mode instability, although only one 
transducer was available for this measurement and the maximum amplitude may not have been measured. Perhaps 
more concrete proof this actually was transverse mode instability is provided by a test of a similar injector but in a 
calorimeter instead of ablative chamber. 

High-frequency instability occurred about 1.6 seconds after initiation of LCH 4 on this test. Figure 38 shows the 
transition from chug to high-frequency instability. The dynamic nature of this transition, starting from the peak of 
the chug oscillation cycle, closely resembles Fig. 34(b). A spectrogram of this instability is shown in Fig. 39. 
Unfortunately, this transducer failed shortly after the time shown. The heat flux in the combustion chamber in this 
test, measured with water-cooled calorimeter spools, increased dramatically within 1.5 seconds after the initiation of 
the instability, as shown in Figs. 40 and 41. The test was terminated just after the time of the latest heat flux average 
shown on Fig. 40 due to high calorimeter outlet water temperature. 

An 8- to 12-times increase in heat flux at the injector face as indicated by Fig. 40, along with nearly a doubling 
of heat flux in the barrel and at the throat, is certainly the result of high-frequency combustion instability at a 
transverse acoustic mode of the combustion chamber. The only possible mode for this frequency is the IT, as 
predicted at a frequency about 90% of that calculated with equilibrium sound speed. Instabilities due to coupling to 
injector features such as oxidizer post modes have not been observed to develop such thermal distress. 

This instability has different characteristics than the short duration spontaneous instability shown in Fig. 32 for 
the 28-element injector - no apparent source frequency such as the 2540 Hz signal is apparent in Figs. 37 or 39. For 
this and the other two high-frequency instabilities, the Hewitt Stability Correlation provides some guidance. Hewitt 
d/V are plotted versus chamber diameter in Figs. 42 and 43 for oxidizer and fuel, respectively, where d for oxidizer 
is film thickness and d for fuel is gap width. 

Several key aspects can be diagnosed from Figs. 42 and 43. Foremost, comparison of the oxidizer d/V using 
oxidizer swirl flow film thickness as d predicts three 40-element injector tests conducted with smaller oxidizer 
metering features that went unstable to be unstable at the IT mode of the undamped 5.66” diameter combustion 
chamber. Tests with 40-element injectors with larger oxidizer metering features, which were chugging and had rich 
high-frequency content but did not have sustained high-frequency instability, are located in a more stable zone 
between the V 2 -IT line and IT line, as shown in Fig. 42. Injectors located between the %-lT line and the IT line 
most likely will require combustion chamber damping devices for stabilization. 35 No damping devices were 
included in any of the 5.66” diameter chamber tests. 

Second, oxidizer d/V for the 28-element injector tests, also located between the %-lT and IT lines, are likely to 
be stable or marginal in IT in undamped combustion chambers, and have increased levels of high-frequency content 
in the dynamic data, which is also exactly what was observed. Also, while the maximum sustainable frequencies 
from oxidizer d/V shown in Fig. 42 for the 28-element injector do not support instability, the proximity to the IT 
line suggests some influence of the IT mode. This is the proposed reason why 4000-Hz signals in lower chamber 
pressure tests (at lower oxidizer velocities) become 4500-5000-Hz signals at higher chamber pressure tests (at higher 
oxidizer velocities). 
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Figure 38. Transition from chug to high-frequency 
instability during Test P8006A-009. 
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Figure 39. Spectrogram of high-frequency instability 
during Test P8006A-009. 
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Figure 40. Increase in measured combustion chamber 
axial heat flux due to high-frequency instability in 


Test P8006A-009. 


Figure 41. Increase in measured combustion 

chamber circumferential heat flux due to high- 
frequency instability in Test P8006A-009. 




Figure 42. Hewitt Stability Correlation for unlike 
doublets, with oxidizer d/V calculated for the 
swirl coax with d as oxidizer film thickness. 


Figure 43. Hewitt Stability Correlation for unlike 
doublets, with fuel d/V calculated for the swirl 
coax with d as fuel gap width. 
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Finally, comparison of the fuel d/V using fuel gap width as d as shown in Fig. 44 predicts the 40-element injector 
tests to be stable, contrary to the evidence. Thus, from a d/V perspective, the LCH 4 circuit is not the driving source 
for instability in these tests. For the 28-element injector, Fig. 43 shows the fuel system is less stable, so in tests of 
that injector the fuel may be generating high-frequency content. 

IV. Summary and Conclusions 

As part of an effort funded by the Propulsion and Cryogenic Advanced Development (PC AD) project, under the 
Exploration Technology Development Program, the NASA Marshall Space Flight Center has been evaluating 
combustion performance and stability modeling capability for oxygen and methane propellants. The purpose of 
conducting these analyses was to evaluate the capability of currently available methods to predict combustion, 
performance and stability characteristics for future development of liquid rocket engine combustion devices with 
methane propellant. 

Combustion performance analyses have been conducted on several rocket engine thrust chamber configurations 
using liquid oxygen and gaseous and liquid methane propellants where characteristic velocity performance could be 
evaluated. For several of these cases, axial heat flux profiles were available as well, which provided an additional 
anchoring technique to the combustion analyses. The results of performance analyses of the selected L0 2 /CH 4 test 
cases with the ROCCID and CICM codes provided mixed results. The ROCCID code provided acceptable 
predictions for single-element shear coaxial element injectors with liquid oxygen and gaseous methane propellants 
at higher chamber pressures (greater than about 600 psia), using the Burick correlation for E m and a 30% increase in 
the DROPMIX droplet size. However, predictions for throttled conditions (low chamber pressures -300 psia), given 
anchoring at higher chamber pressures, were unacceptable. Operating conditions at the low chamber pressures were 
within the range of a mixing correlation for shear coaxial elements, but use of the mixing correlation significantly 
over-predicted performance. Additional mixing, atomization and hot-fire test data are necessary to examine the 
discrepancies at these conditions. 

The CICM code also provided reasonable results for single-element shear coaxial element injectors with liquid 
oxygen and gaseous methane propellants at higher chamber pressures when new atomization anchoring parameters 
were used. The manual-suggested inputs for cup and chamber atomization parameters, which were developed from 
oxygen and hydrogen propellant engines, were not valid with oxygen and methane propellants. Development of 
new atomization parameters with recommended methods was not successful since cup pressure drops in the 
L0 2 /gCH 4 test cases were negligible. Consequently, parametric analyses were conducted to determine a set of 
parameters that could best fit the test data. Two of the four atomization parameters were found to be influential. 
Calculations with these parameters provided reasonable comparisons to both performance and axial heat flux profile 
data, but like with the ROCCID cases, only at higher chamber pressures (greater than about 600 psia). At the low 
chamber pressures for throttled conditions (300 psia) where ROCCID over-predicted measured performance, 
analyses with CICM greatly under-predicted the measured performance. 

With swirl coaxial elements with liquid oxygen and gaseous methane propellants, significant alterations of 
atomization parameters calculated within ROCCID and E m from available correlations were required to match 
measured performance and energy release rates suggested by axial heat flux profile measurements. Significant 
improvements in mixing and atomization predictions are required for swirl coaxial element injectors. 

For cases with liquid methane with multi-element LOL impinging element injectors, ROCCID provided 
acceptable predictions using the Nurick correlation for E m and the DROPMIX atomization parameters as calculated 
in ROCCID. With swirl coaxial elements, significant changes were again required for the calculations of 
atomization parameters within ROCCID. The atomization calculation for both the swirled liquid oxidizer flow and 
the annular liquid fuel flow required extensive modifications to fit to the measured data, and failed outright for 
elements with large fuel annulus widths. However, use of quasi-anchored parameters predicted the vaporization- 
limited performance that likely occurred with several of the cases. Significant effort is required to improve the 
capability to make such predictions. 

Combustion stability characteristics of several multi-element swirl coaxial element injectors with L0 2 /gCH 4 
propellants were examined and correlated with simple models. These injectors were designed, fabricated and tested 
by the NASA MSFC with funding from PC AD. Low frequency combustion instability (chug) occurred with most 
injector patterns, with frequencies ranging from 150 to 250 Hz. While several tests with a 28-element injector had 
low peak-to-peak chug oscillation amplitudes (0.7%-4% of mean chamber pressure), no test appeared to be 
completely free of chug influence, even at chamber pressures exceeding 360 psia and mixture ratios as low as 2.6. 
Modeling with the methods of Wenzel and Szuch showed the influence of chamber pressure, mixture ratio and 
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LCH 4 injector pressure drop on neutral chug stability. Increasing LCH 4 injector pressure drop at the injector 
element exit was predicted to improve chug stability margin for the 28-element injector, although stabilization at 
desired operating conditions will require increased injector pressure drops on both circuits. Tests with a 40-element 
injector with low LCH 4 injector pressure drop at the injector element exit exhibited similar chug characteristics. A 
later revision of the 40-element injector with increased LCH 4 injector pressure drop at the injector element exit as 
well as increased L0 2 injector pressure drop appeared to improve chug stability, although reliable dynamic data 
were not available for these tests. 18 

All tests with the 28-element swirl coax element injector and the handful of tests with the initial 40-element swirl 
coax element injector patterns that included a high-frequency pressure sensor showed a rich high-frequency dynamic 
environment in the undamped combustion chambers. The chug oscillations for these injectors were contaminated by 
high-amplitude, high-frequency ringing at the maximum pressure of the chug oscillation cycles, with some bursts 
lasting several tens of milliseconds, suggesting such injector patterns might be unstable in dynamic stability rating 
tests in these undamped combustion chambers. 

High-frequency chamber pressure dynamics with the 28-element injector appears to be the result of two 
phenomena - oscillations of the gaseous core in the oxidizer post, and high-frequency instabilities due to unsteady 
atomization of the conical oxidizer film. For the former, simple estimates of “organ pipe” frequencies of the 
gaseous core compared reasonably well to some observed oscillation frequencies. Spontaneous bursts of high- 
frequency activity from such an oxidizer post mode may have excited the IT mode on several tests, exceeding 
allowable damp times for the IT mode. Mechanistic models to predict the responses of such oxidizer post “organ 
pipe” oscillations for a hollow-core swirl coaxial element are in development. Note that while this post resonance 
may have been destabilizing in the current case, the use of a properly tuned gas core resonance as a stabilizing 
feature is also possible. 47 For the atomization-dependent phenomenon, the Hewitt Stability Correlation provided a 
reasonable rationale for low-amplitude excitement of near- 1 T mode frequencies that changed with chamber pressure 
(actually, oxidizer injection velocity). 

High-frequency chamber pressure dynamics with the initial 40-element swirl coax element injector patterns for 
which a high-frequency pressure sensor was available appeared similar to the 28-element injector dynamics, except 
that the 40-element injector was more sensitive to the IT mode. High-frequency combustion instability occurred in 
three tests apparently at the IT mode with initial versions of the 40-element injector, evident by a frequency 
measured with a single high-frequency chamber pressure measurement as well as a dramatic increase in head-end 
chamber heat flux measured in one test with a calorimeter chamber. The occurrences of these instabilities correlated 
closely with predictions from the Hewitt Stability Correlation using oxidizer film thickness as a correlating 
parameter. These results suggest that the 40-element injector pattern (and most likely the 28-element injector 
pattern as well) would require combustion chamber damping devices for high-frequency stabilization. 

V. Recommendations for Future Work 

It is apparent that a priori predictions from the non-CFD codes of combustion, performance and stability that 
were examined with methane (and especially liquid methane) propellant are not in a reliable state to support design 
of new combustion devices with this propellant. The primary errors appear to extend from failure to predict 
atomization characteristics. Some element types are in a better state of readiness than others, e.g., shear coaxial 
elements with gaseous methane and like-on-like doublet elements with liquid methane. Also, in several cases when 
anchored to hot-fire test data, the existing codes can predict combustion performance reasonably well, with the 
exception of deep throttled conditions. For high quality predictions, data for similar geometries are needed to 
anchor the prediction. 

However, improvement of capability to calculate a priori combustion, performance and stability will be 
necessary for the cost-effective development of new combustion devices for liquid rocket engines using methane. 
The following list provides recommendations for the improvement of combustion performance analysis modeling 
techniques with L0 2 /gCH 4 and L0 2 /LCH 4 propellants. 

1) Extend the Burick shear coaxial element Rupe mixing efficiency (E m ) correlation for conditions with 
methane fuel. Use of an element with the higher densities of methane (or any fuel denser than gaseous 
hydrogen) required extrapolation beyond the range of the correlation in all cases except at the very low 
chamber pressures, when it appeared that the mixing efficiency from Burick was significantly 
overpredicted. The mixing efficiencies could be derived experimentally or computationally. 

Calculating the mixing efficiency numerically would be acceptable, given the predictions are validated 
against experimental data already derived. 
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2) Measure and correlate mixing efficiency for liquid/gas swirl coaxial element flows similar to the 
Burick correlation. Since there are little data currently available, this should include an experimental 
program, although conducting numerical experiments is also recommended as in #1. 

3) Extend or improve mixing efficiency correlations for liquid/liquid swirl coaxial element flows. 

4) Improve the atomization model for swirled liquid injected from the center orifice of a conventional 
swirl coaxial element. This improvement is a key need for analyzing swirl coaxial elements for any 
injector configuration that uses this element. 

5) Generate atomization data for swirled liquid injected from the center orifice of a swirl coaxial element, 
including breakup length data and drop size data. These data will be critical to validating the 
atomization model described in #4. 

6) Improve the atomization model for liquid injected from the outer annulus of a coaxial element. Both 
axial and swirl injection should be considered. This model will make analysis of the liquid/liquid 
coaxial element feasible. Predictions from currently available models of this flow are unreliable. 

7) Generate atomization data for liquid injected from the outer annulus of a coaxial element, including 
breakup length data and drop size data. These data will be critical to validating the atomization model 
described in #6. 

8) Ensure that the mixing and atomization data are generated over a wide range to accommodate deep 
throttling conditions. 

9) When conducting hot- fire test programs, extend the tested range of conditions, especially to lower 
chamber pressures, to generate a wide range of parametric data for improving model capability. 

For stability modeling, further development is required to calculate organ pipe resonance frequencies and their 
responses in the oxidizer posts of swirl coax elements, along with further elucidation of the influence of d/V - type 
instabilities. In addition, with liquid methane, stabilization devices in the combustion chamber will almost certainly 
be required when such conventional elements and designs are used, so further development of modeling the 
influences of damping devices will be necessary. 

This paper does not include results from calculations with combustion computational fluid dynamic (CFD) 
models. The authors certainly recognize that future capability with combustion CFD models will eventually 
overtake the one-dimensional models. However, the mixing and atomization model developments recommended 
above will be useful (or indeed, critical) in development and validation of combustion CFD models as well. 

In the interim before the general use of CFD for complete modeling of performance and stability, CFD models 
can be used in other ways. Use of CFD to calculate parameters of importance to one-dimensional (non-CFD) 
modeling has been under development at the NASA MSFC for some time. Recently, calculations of E m of shear 
coaxial element flows with F0 2 /gCH 4 propellants have been conducted using Foci-CHEM. Such calculations can 
provide significant improvement to one-dimensional model inputs and possible preclude the necessity of measuring 
such parameters in difficult-to-run experiments. 
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Motivation - Technology Development of 
L0 2 /CH 4 LRE Combustion Devices 


Marshall Space Flight Center 


• NASA considering use of oxygen and methane 
propellants for future in-space missions 

- Improves performance over storable propellants 

- Provides compatible oxidizer (oxygen) with life support systems 

- Provides possibility for in-situ propellant production 

• Oxygen-methane propellants are not used on any 
current flight-qualified space propulsion system 

• Significant increase in Technology Readiness Level 
(TRL) of thrust chamber technology is required 

- Capability to analyze thrust chamber phenomena 

• Combustion, performance, heat transfer, stability, ignition 
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Objectives 


of NASA Combustion Devices 
Analysis Plans 


Marshall Space Flight Center 


• Long Term: 

- Improve and validate capability to analyze multi-element thrust 
chamber combustion with oxygen and methane propellants 

• Performance 

• Combustion Stability 

• Heat transfer 

• Transients 

• Near Term: 

- Evaluate current capability to analyze combustion, performance 
heat transfer, and stability 

• L0 2 /gCH 4 

• L0 2 /LCH 4 

• Single element injectors 

• Multi-element injectors 
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Evaluate Capability to Predict Combustion 
Performance, and Heat Transfer 


Marshall Space Flight Center 


• Task 1 - Collect information for several test cases 

- Penn State/CUIP LOX/gCH4 single element coax (2006) 

• Shear coax - recessed oxidizer post 

• Shear coax - flush oxidizer post 

• Swirl coax - flush oxidizer post 

- Rocketdyne LOX/gCH4 82-element shear coax (1989) 

- NASA MSFC LOX/LCH4 40-element like-on-like paired doublet (2007) 

- NASA MSFC LOX/LCH4 28-element swirl coax (2008) 

- NASA MSFC LOX/LCH4 40-element swirl coax (2008, 2009) 

• Task 2 - Analyze all cases with heritage 1-D codes: 

- Rocket Combustor Interaction Design and Analysis (ROCCID) 

- Coaxial Injector Combustion Model (CICM) 

• Task 3 - Analyze progression of cases with combustion CFD 

- Loci-CHEM 

- Not discussed in this paper 
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Marshall Space Flight Center 


Penn State University/CUIP L0X/gCH4 
Single-Element Coax Injectors 
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Shear Coax Injector 
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Locke, J.M., Pal, S., and Woodward, R.D., “Chamber Wall Heat Flux Measurements for a 
JACOBS LOX/CH4 Uni-element Rocket,” AIAA 2007-5547, July 2007. 
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PSU CUIP L02/gCH4 Single-Element Tests 
Measured Relative Performance 

Marshall Space Flight Center 



Injector-End Chamber Pressure, psia 
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One-dimensional Combustion Models 


Marshall Space Flight Center 


• Rocket Combustor Interactive Design and Analysis Program (ROCCID) 

- Developed at Aerojet for NASA GRC in 1991 

- Commercial version currently available from Sierra Engineering 

- Capable of evaluating many element types 

• like-on-like doublets, shear coax, swirl coax (gas and liquid fuels) 

- Steady-State Combustion Iteration (SSCI) 

• Propellant atomization calculated from mechanistic models but atomization parameters can be altered 
by user 

• Vaporization calculated with Priem & Heidmann Generalized Length correlation 

• Mixing calculated using two-zone four-streamtube model 

• Coaxial Injection Combustion Model (CICM) 

- Developed at Rocketdyne for NASA MSFC in 1974 

- JANNAF-standard computer program for evaluation of liquid rocket engine performance 

- Capable of analyzing shear coaxial elements only 

- Steady-state combustion routine 

• Calculates rate of atomization from shear force between jet and surrounding gas 

- Calculates an axially varying mean droplet size 

- Does not require experimental dropsize correlations but does require user-defined parameters 

• Specially developed droplet heating and vaporization model 

• Provides significant mechanistic calculations of combustion in the recessed oxidizer post region - the 
“cup” 

• Intra-element mixing (Em) defined by user - run and average multiple runs at different mixture ratios 

- Then interface with a Distributed Energy Release (DER) or Streamtube Combustion (STC) 

ESTS Group code 
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Methane Shear Coax Elements Out of Range 


Marshall Space Flight Center 


of Burick Liquid/Gas E m Mixing Correlation 
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2008 ROCCID Calculation - 
Recessed-Post Single-Element Shear Coax 
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Vapor O/F or Heat Flux, BTU/in A 2-sec 
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ROCCID Calculation of PSU Shear Coax at 
Pc - 300 psia Has Uncertainties 
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Vapor O/F or Heat Flux, BTU/in A 2-sec 



With New L0 2 /gCH 4 Atomization Parameters, CICM OK 

at High Pc But Not Good at Low Pc 
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Could not use recommended methods - 
no measureable cup AP with methane 
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NASA MSFC Conducted Combustor Testing 
with LOo/LChL Propellants For Several Years 
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1 . Elam, S., and Sprader, C., “Test Summary Report for Test Program 
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LCH4/L0X ROCCID Calculation with Like-on-Like 
Paired Doublet Multi-element Injector is Good 
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28 Element ROCCID Calculation of 
Performance Change with L’ is Reasonable 
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User-Inputs for Atomization of 28 Element Liquid/Liquid 
Swirl Coax with Small Fuel Annulus Gaps 
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Chug Waveform with 28-Element Injector 
W Contained Lots of High-Frequency Content 
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Oscillations at 2500 Hz May Be Exciting IT Mode 

Marshall Space Flight Center 

• Closer look shows extended high frequency oscillation 
periods 

- some more than “allowable” for acoustic modes 
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Damage at Oxidizer Swirl Inlets Was 

Observed 
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High-Frequency Instability Initiated 
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High-Frequency Instability Occurred Again 

In Calorimeter Chamber 
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• 4490 Hz is 89% of equilibrium f 1T 

• Factor of 1 1 increase in head-end heat 
flux 

• Factor of 1 .8 increase in throat heat flux 

...very likely the IT mode 
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4000 Hz and 5000 Hz Oscillations may be 

from Hewitt “d’VV 


Marshall Space Flight Center 


• Hewitt Stability Correlation for “unlike doublets” 

- Use for 
propellant 
impingent 

- “Twice as 
as like-on-li 
doublets 

• For fuel “d” (with liquid methane) use annulus 
gap width 

• For oxidizer “d” use swirl film thickness (from 
COBRA preburner presented in 2004 JANNAF) 
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Hewitt Stability Correlation Predicts IT Mode 
*Sr Instability with Oxidizer “d”N 
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Summary and Conclusions 


Marshall Space Flight Center 


• a priori performance predictions with 1-D codes with methane (and 
especially liquid methane) are not adequately reliable to support design of 
new combustion devices 

- Significant adjustments to atomization inputs for most cases 

- Predictions of liquid/gas shear coax ok at higher Pc but not at low Pc 

- Predictions with liquid/Liquid swirl coax show vaporization-limited combustion although 
with extensive user intervention 

• Low-frequency combustion instability (chug ) 

- Analysis models are adequate 

- Need further evaluation for why inlets broke on L0 2 /LCH 4 swirl coax injector 

• High-frequency combustion instability 

- Analysis of L0 2 /gCH 4 shear coax element injectors focus on oxidizer post resonances 

- Analysis models for L0 2 /LCH 4 swirl coax injectors need work 

• Oxidizer “gas core” resonances 

- Gas is combustion products 

- Tune it out or tune it in (“IT resonator” - Cha et.al, AIAA 2005-3748) 

• IT mode instability driven by Hewitt “d”A/ of oxidizer film thickness 

- Liquid oxygen side, not the liquid methane side 

- Previously observed during COBRA subscale preburner testing 
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Recommendations 


Marshall Space Flight Center 


1 ) Extend the Burick shear coaxial element Rupe mixing correlation 
for methane (higher density) fuel. 

2) Improve the atomization model for liquid methane injected from the 
outer annulus of a coaxial element. 

3) Improve the atomization model for swirled liquid oxygen injected 
from the center orifice of a swirl coaxial element. 

4) Generate atomization data (cold flow and hot-fire) for the 
liquid/liquid swirl coaxial element, including breakup length data 
and drop size data. 

5) Develop new methods to calculate the atomization parameters for 
the CICM model. 

6) Extend the tested range of conditions during hot-fire test programs, 
especially to lower chamber pressures, to generate a wider range 
of parametric data for validating model capability. 
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